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Abstract. Grobner Bases |Buc70] and Cylindrical Algebraic Decomp- 
osition [C0I75ICMMXYO9] are generally thought of as two, rather differ- 
ent, methods of looking at systems of equations and, in the case of Cylin- 
drical Algebraic Decomposition, inequalities. However, even for a mixed 
system of equalities and inequalities, it is possible to apply Grobner bases 
to the (conjoined) equalities before invoking CAD. We see that this is, 
quite often but not always, a beneficial preconditioning of the CAD prob- 
lem. 

It is also possible to precondition the (conjoined) inequalities with respect 
to the equalities, and this can also be useful in many cases. 



The examples used in this paper are available in jWill2| . This work was partially 
supported by the U.K.'s EPSRC under grant number EP/J003247/1. 

1 Introduction 

Solving systems of equations, or equations and inequations (7^) /inequalities 
(>, <) is an old subject. Deciding the truth of, or more generally eliminating 
quantifiers from, quantified Boolean combinations of such statements, is more 
recent [TarSl, . We can distinguish many families of methods, even if we restrict 
attention to the real numbers, or possibly the complex numbers. 

=G The method of Grobner bases. Here the input is a set S — {si, . . . , Sk} 
of polynomials in some polynomial ring k[xi, . . . , Xn] equipped with a total 
ordeilj -< on the monomials, and the output is a set G = {pi, . . . ,pi} which 
is equivalent, in the sense that it generates the same ideal, i.e. (G) — (S), 
and is simpler, or "surprise-free" , in that the leading monomial with respect 
to ~< (denoted Im^) behaviour is explicit, (Im^(G)) = (lm^((G))). Then the 
solutions of G are those of S, i.e. 

{x : pi(x) = A P2(x) - A • • • Ap,(x) = 0} . (1) 



^ We have concentrated on purely lexicographical orders, since these seem to be the 
most useful to us. 



=A The method of triangular decomposition via regular chains |ALMM99IMM05"] . 
Here the output is a set of regular chains of polynomials 



{(Pl,l'-Pl,2, ■ • {P2.1,P2.2, ■■■),■■ •}, 



(2) 



and the solution is the union of the set of regular zeros of these regular 



<Coi The method of Cylindrical (semi-)Algebraic Decomposition for real closed 
fields, computed via repeated projection to and repeated lifting [Col75[ 
and many improvements]. 

=Coi The previous case restricted to equality. 

<AR The method of Cylindrical (semi-)Algebraic Decomposition for real closed 
fields via triangular decomposition |CMMXY09| . 

The method of Cylindrical Decomposition over the complexes via triangu- 
lar decomposition, which was introduced in |CMMXY09| as a stepping-stone 
to the previous method, but which probably has independent interest. 

/9cH Quantifier Elimination by partial (i.e. taking account of the Boolean struc- 
ture and quantifier structure) Cylindrical Algebraic Decomposition |CH91) . 

Others such as Weispfenning's Virtual Term Substitution |Bro05[ is a readable 
introduction], or Tarski's original method jTarSl^. 

Conversely instead of asking for solutions x to 3x/i (x) > A • • • , we may use 
a Positivstellensatz to show that no such x exist, as in |PQR09| . We do not 
discuss this direction further here. 

It should be noted that both <coi and <aii/ y^AC (but not /3ch) have the draw- 
back that the Cylindrical Algebraic Decomposition produces decompositions for, 
not only the question posed, e.g. yy3zp{x, y, z) — OAq{x, y, z) = OAr{x, y, z) > 0, 
but also all other questions involving the same polynomials, provided the quan- 
tifiers are over variables in the same order, e.g. 3yVzp{x, y,z) < {q{x, y,z > 
Ar{x,y,z) = 0). 

This paper asks the question: "can these methods usefully be combined?" 
The combinations we are thinking about are those of conjunction: Can the fact 
that B is in the context of ai = A • • • A = A i? be used to simplify Bl In 
particular, we look at the use of Grobner base methods to simplify the equalities 
in the conjunction and to simplify the inequalities in the light of the equalities. 

Technical Note: all computations (=Gi <aii and j^ac) were performed in 
Maple 16/3 on a 3.1GHz Intel processor, except for the <coi, =Coi and ^ch ones, 
which were performed on a 2.83GHz Intel processor with QEPCAD B version 
1.65 |Bro03) . Times for a hybrid calculation, e.g. =g/<Co1, are either quoted as 
the total time or a decomposition a + b — c where a is the time (in milliseconds) 
for —a, b for <coi, and c is the sum. We have run QEPCAD in three modes: 



chains: 





1. on the problem as given in |BH91) . implementing ^cW, 



2. as above but with the full-cad option to ignore the Boolean structure of 
the expression; 

3. with no quantifiers stated, and the full-cad option, implementing <coi- 

2 Examples in this paper 

2.1 |BH9l| 

This paper has a variety of examples for ^ca, all of a form to which =g is 
applicable. 

2.2 |CMMXY09l 

This paper has a variety of examples for <ziR. We chose some of those to which 
=G is applicable. 

2.3 Two Spheres and A Cylinder 

Let the following be spheres in R^: 

Si: (a;-l)2+y2 + ^2_3. 

^2 : {x + 1)2 + y2 + ^2 _ 3. 

53: {x-lf + {y-^y + z^~i; 

S,: (a; + l)2 + (y+|)V(z + |)2-3. 

Denote the infinite cylinder centred on the z-axis with radius 1 by C, so that 
the equation defining the cylinder is: 

C: 

Now we investigate intersecting pairs of spheres (roughly increasing in CAD 
'difficulty') under conditions based on the cylinder. We assume the spheres' 
equation will always be required to equal but make no assumptions on the 
condition on the cylinder. That is, we wish to solve the problem: 

5, -0 A -0 A C*0 *£ {=,^, <,>,<,>}, i = 1,2,3. (4) 

We use the underlying variable ordering (2, y, x). 

^ This is the QEPCAD notation, meaning that we will project from (z, j/, a;)-space 
to (z,y)-space to (z)-space. We therefore end up with polynomials in z alone, 
so this is equivalent to a purely lexicographical Grobner base with z < y ^ 
X, i.e. plex( [z,y,x] ) in Maple: =gc is used to indicate Grobner bases with 
this (compatible) ordering. The CAD package in Maple [CMMXY09| requires 
PolynomialRingC [x,y,z] ) to achieve the same effect as QEPCAD's {z,y,x). —qr 
denotes the reverse plex order. 



3 Prior Art 



Needless to say, we are not the first to have had this idea. 
3.1 Buchberger Hong 

|BH91) considers the case of =g f [BGK85] re-implemented in C) applied to <coi 
(an early version of |CH91j re- implemented in C), i.e., rather than computing 
a CAD for the zeros of a system of equations E (i.e. ei = A 62 = A • • • ) 
and inequalities F, compute it for G, a (purely lexicographical) Grobner base 
for i?, and F. They generally found a very substantial speed-up in the total 
computation time, e.g. "Solotareff A'H 

3x3y ix^ — 2x — a = — — ax — 2b + a — 2 = (5) 
3?/^ - 2y - a = - 2/^ - ay - 26 + a - 2 = A (6) 
4a e [1, 7] A 46 G [-3, 3] A a; £ [-1, 0] A y £ [0, 1] (7) 

(with the variable ordering {b,a,x,y)) took them 11500 ms for ^Sch, but 717 
for =G, and 117 for ^ch applied to the result, a total of 834 ms, or a 13- 
fold speed-up. "Solotareff B" is the same problem but with (a, b, x, y) as the 
variable ordering, and here the /)ch time was again greatly reduced, but the =g 
time was excessive. Of course, there have been substantial improvements in the 
implementation of all these algorithms since [BH91| was published, and Table [2] 
shows that the =g time is now less than 1/3 of the /3ch time. We choose rather 
to focus on the number of cells generated, which is closely connected to the =coi 
time, and also affects the time taken to make use of the output. The cell counts 
are shown in Table [TJ 



Table 1. Cell counts for Solotareff 





Ordering A 

<Col =g/<Co1 


Ordering B 

<Col =g/ <Co1 


dSHZl Partial 
Full 


153 63 
349 625 


375 41 
1063 237 


(EHSI Partial 
Full 


29 15 
29 33 


97 17 
97 17 



More reruns of |BH91) are given in Table [5J We see that, with today's tech- 
nology, the conclusion of |BH91| . viz. that =g generally improves /9ch for the 
class of problems to which it is applicable, is still generally valid, but the details 
differ: notably the Grobner base time is generally insignificant today. 



^ There are various problems labelled "Solotareff": for a description of this class see 
[Will2j and the links therein. 



Table 2. |BH91j with today's technology 









=G 


/^CH 




^cn/full-cad 


=G/;acH/full-cad 




Time Cells 


Time 


Cells 


Time 


Cells 


Time Cells 


I A 


190 


503 


22+72= 


94 


23 


188 


503 


22+73= 


95 51 


I B 


199 


369 


21+74= 


95 


17 


191 


369 


21+75= 


96 33 


R A 


85 


1 


24+73= 


97 


1 


86 


1 


24+71= 


95 1 


R B 


129 


1 


24+72= 


96 


1 


125 


1 


24+72= 


96 1 


E A 


297 


621 


25+134= 


159 


621 


576 


11139 


25+394= 


419 11139 


E B 


Error 


? 


50+?= 


Error 


? 


Error 


? 


50+?= 


Error ? 


S A 


89 


153 


22+72= 


94 


63 


199 


349 


22+185= 


207 625 


S B 


113 


375 


23+75= 


98 


41 


228 


1063 


23+180= 


203 237 


C A 


133 


19 


42+?= 


Error 


? 


235 


19 


42+?= 


Error ? 


C B 


Error 


? 


132+?= 


Error 


? 


Error 


? 


132+?= 


Error ? 



Table 3. IBH91! Examples for full CADs 







Col 


=g/=Co1 






=g/</1R 




Time 


Cells 


Time 


Cells 


Time Cells 


Time Cells 


I A 


236 


3723 


22+77= 99 


273 


29426 3763 


2470 273 


I B 


212 


3001 


21+76= 97 


189 


36262 2795 


1482 189 


R A 


150 


2101 


24+86= 110 


105 


17355 1267 


570 165 


R B 


21091 


7119 


24+80= 104 


141 


356670 7119 


470 141 


E A* 


7390 114541 


25+3189= 3214 53559 


262623 28557 


62496 14439 


E B* 


Error 


? 


50+?= Error 


? 


> 1000s ? 


> 1000s ? 


S A* 


115 


1751 


22+82= 104 


297 


16014 1751 


2025 297 


S B* 


253 


6091 


23+82= 105 


243 


43439 6091 


1647 243 


C A* 


820 


8387 


42+?= Error 


? 


216028 7895 


> 1000s ? 


C B* 


Error 


? 


132+?= Error 


7 


> 1000s ? 


> 1000s ? 



indicates that the linear inequalities have been omitted in this version. 



There is one point which is not exphcit in jBH91| . As the computation of 
Grobner bases in one variable is just equivalent to Euclid's algorithm, i.e. Gaus- 
sian elimination in Sylvester's matrix, Grobner base computations which are not 
genuinely multi-variate do not affect the set of resultants etc. generated in <coi, 
and hence are of limited use in the projection phase. They might still reduce the 
work done in the lifting phase, of course. 

Table [3] re-runs the examples of |BH91j , but asking for complete cylindrical 
algebraic decompositions, and hence we can compare <coi with <aii legiti- 
mately. Given that the algorithms are fundamentally different, the similarities 
in cell counts are striking. The differences in cell counts (where present) re- 
flect differences in the cylindrical algebraic decompositions for the same input 
problem. 



3.2 Phisanbut 



Phisanbut [Phillj , considering branch cuts in the complex plane, observed that 
g = A f > could be reduced to g = A prem(/, g) > under suitable 
conditions, where prem denotes the pseudo-remainder operation. More precisely, 
if / and g are regarded as polynomials in the main variable x, of degrees d 
and e respectively, then prem(/, g) — rem{c'^~^^^ f , g) , where c is the leading 
coefHcients of g. When g — and c > 0, or when d — e + 1 is even, prem(/, t;) 
has the same sign as /. Unfortunately c might have variable sign, and d — e + 1 
might be odd, so define pprccond(/, g) = rcm(c'-''~'^"''^^ f,9), where n* is n if n is 
even and n + 1 if n is odd. Maple also defines sprem(/, 5) = rem(c™/, 5), where 
m is the smallest integer such that the division is exact, and by analogy we have 
sprecond(/, 5) — rem(c™ f,g)- Note that sprecond(/, 5) = pprecond(/, g) or a 
strict divisor of it, i.e. sprecond is never worse. She generally, but not always, 
saw [Phil 11 Tables 8.13, 8.14] a significant decrease in the number of cells, and 
the time taken to compute sprecond was minimal. 

4 Further developments 

4.1 —a with </iR, 

It would seem natural to apply —q to </iR, as |BH91) did to ^ch- The results are 
in Table 131 and show a speed-up in all instances except the Collision problems. 
We also note the substantial speed advantage enjoyed by <Coi, and this is a 
subject for further study. 

4.2 —G with ^/ic 

We can also mix —g with =/=ac, and these results are shown in Table |4l which 
also compares ^ac with <ziR. <ziR involves doing ^ac first, and then running 
the MakeSemiAlgebraic algorithm from [CMMXY09] . For these examples, the 
MakeSemiAlgebraic step is the most expensive initially, but often not after we 
apply =G- 

4.3 =G with inequalities in <^r 

Having reduced the equalities to a Grobner base G, it is now possible to reduce 
the inequalities by G, since adding/subtracting a multiple of an element of G is 
adding/subtracting 0. We can reduce with respect to the main variable, denoted 
=g/^xj ■^ith respect to secondary variables, denoted =G/^y , or with respect 

to all variables (Maple's NormalForm), denoted =0/-^ ■ If we compare tables 
[6] and [7] we see that the number of cells produced is the same across the two 
methods. 



Table 4. Timings for |BH91) Examples: <z1r//z\c 







<ZiR 












Time 


Time 


Ratio 


Time 


Time 


Ratio 


Intersection A 


5691 


29426 


4.17 


1168 


2470 


1.11 


Intersection B 


5584 


36262 


5.49 


886 


1482 


0.67 


Random A 


4614 


17355 


2.76 


310 


570 


0.84 


Random B 


67343 


356670 


4.30 


318 


470 


0.48 


Ellipse A* 


85425 


262623 


2.07 


27916 


62496 


1.24 


Ellipse B* 


441245 


> 1000s 




> 1000s 


> 1000s 




Solotareff A* 


6666 


16014 


1.40 


1760 


2025 


0.15 


Solotareff B* 


9536 


43439 


3.56 


1404 


1647 


0.17 


Collision A* 


41085 


216028 


4.26 


> 1000s 


> 1000s 




Collision B* 


> IQQQs 


> 1000s 




> 1000s 


> 1000s 





"Ratio" = (<ziR— 7^zic)//zic, i.e. the relative cost of MakeSemiAlgebraic. 



Table 5. Examples from [CMMXY09] 





<AR. 






AR. 




Ratio 




Time Cells 


Time 




Cells 


Time Cells 


Cyclic-3 


3136 


381 


20 + 245 = 


265 


21 


11.83 18.14 


Cyclic-4 


> 1000s 


? 


64 + 5813 = 


5877 


621 


? ? 


2 


2249 


895 


22 + 1845 = 


1867 


579 


1.20 1.55 


4 


3225 


421 


24 + 19738 = 


19762 


1481 


0.16 0.28 


6 


363 


41 


20 + 918 = 


938 


89 


0.39 0.46 


7 


3667 


895 


26 + 6537 = 


6563 


1211 


0.56 0.74 


8 


3216 


365 


21 + 174 = 


195 


51 


16.49 7.16 


13 


14342 


4949 


18 + 220 = 


238 


81 


60.26 61.10 


14 


334860 27551 


21 + 971 = 


992 


423 


337.56 65.13 



Table 6. Spheres and Cylinders: <ziR 





</lR 

Time Cells 


=g/</1R 

Time Cells 


~g/ / <AU. 

Time Cells 


= g/— /<ZiR 

Time Cells 


= G/A''/<ziR 

Time Cells 


Si,S2,C 
82,83,0 
83, 84, c 


9830 1073 
187048 12097 
247458 11957 


1057 267 
5880 1299 
8164 1359 


394 91 
3171 627 
9177 1123 


528 183 
2149 517 
5476 881 


298 99 
506 213 
590 213 



Table 7. Spheres and Cylinders: <coi 





<Col 

Time Cells 


= g/<Co1 

Time Cells 


1*0. 

=g/^ /<Co1 

Time Cells 


81,82,0 
82, 83, c 
83, 84, c 


30 1073 
763 12097 
1760 11957 


23 + 8 = 31 267 

27 + 36 = 63 1299 

28 + 37 = 65 1359 


24 + 4 = 28 99 

28 + 13 = 41 213 

29 + 14 = 43 213 



5 Choice of Method 



Suppose we are given a problem, which we may formulate as 

quantified variables ei=OA---AeA;=OA . . . , /;), (8) 

where B is a Boolean combination of conditions = 0, 7^ 0, < etc. on some 
polynomials fj, then we may be able, by applying Grobner techniques to the 
Cj , producing , and then reducing the fj , to produce various alternative 
formulations 

quantified variables e^'^ = A ■ • • A e^f,, = A B{f["^ , . . . , //'^ ) , 0''>) 

and each of these may have several variable orderings compatible with the con- 
straints implied by the quantification (if any). Which should we choose? Of 
course, in the presence of arbitrary parallelism, we can start them all, and ac- 
cept the first to finish, but we may wish to be less extravagant. 

In the contexts of ^ch (strictly speaking, the REDLOG implementation), 
and where the only choice was in the variable order, this question was considered 
by [DSS04| . Retrospectively, there are two measures for the difhculty of a CAD 
computation: the time taken and the number of cells produced. For a given ^cu 
problem, they observed that two are usually correlated for different formulations, 
and we observe the same here for <aii — see our tables. However, we would 
like a measure that could be calculated in advance, rather than retrospectively. 

The processes of |Col75ICH91j starts with a set An of polynomials in n 
(ordered) variables xi, . . . ,Xn, and 

f . repeatedly project Ai into Ai^i in one fewer variable, until Ai has only one 
variable, 

* (denote the set {An, . . . , Ai} by A{xi, . . . , x„)) 

2. isolate the roots of these polynomials to get a decomposition of R^, 

3. repeatedly lift the decomposition until we get a (partial for [CH91j ) cylin- 
drical algebraic decomposition of R". 

The third step is, both theoretically and practically, by far the most expensive. 
Hence the question arises: what can we measure at the end of step 1, i.e. depend- 
ing on A only, which is well-correlated with the final cost? Three things come to 
mind. 

ca.Td{A{xi,. . . ,Xn)) =Er=il^»l- 

td{A(xi, . . . ,Xn)) = X]"=i jEAi '^'^(Pi.j) where td denotes total degree. 
sotd(A(xi, . . . , x„)) = E"^, Ep,,,eA, Emonomials m of p., *^(™)- 

|DSS04] discard td, observing that td and sotd are highly correlated and sotd 
"has the advantage of favouring sparse polynomials" . They then observe that 
sotd(A(xi, . . . , Xn)) is significantly more correlated with the retrospective mea- 
sures for any given problem than card. This gives a first algorithm for deciding 
how to project: for all admissible (i.e. compatible with the quantifier structure, if 



any) permutations vr of (xi, . . . , a;„), compute ^(^^^(i), . . . , x^(„)), and choose the 
one with the least sotd value. The drawback of this is that it requires potentially 
(n — l)n\ projection operations. They show that (at least on their examples) this 
always produces a good projection order, and frequently the optimal. 

Table 8. Spheres and Cylinders: <ziR — choice of orderings 







<ziR 


=g/<ar 




/<AR 






Time Cells 


Time Cells 


Time 


Cells 


Si, 5*2, C 


C 


8654 1073 


905 267 


270 


99 




R 




902 267 


453 


183 


S2, S3, C 


C 


189202 12097 


5911 1299 


499 


213 




R 




18941 2639 


5307 


859 


S3, S4, C 


C 


248340 11957 


8159 1359 


580 


213 




R 




160171 9091 


196714 


11203 



They therefore propose a greedy algorithm, where for all permissible choices 
of the first variable to be projected, we compute sotd(v4„„i), and choose the 
variable which gives the least value. Having fixed this as the first variable to 
project, for all permissible choices of the second variable to be projected, we 
compute sotd(^„„2)i and choose the variable which gives the least value, and 
so on. Hence, assuming all projection orders are possible, the number of pro- 
jections done is n + (n — 1) + ■ • • = O(n^) rather than n\. It is currently an open 
question whether the cost of projections behaves similarly. 

We proposed taking this idea still further, and suggested that, for several 
different formulations An, Bn, ■ ■ ■ of a problem, we should compute sotd(A„), 
sotd(i?„), . . . and take the formulation that yields the lowest sotd. We observed, 
however, that neither td nor sotd are good predictors in Table [TTl despite seem- 
ing useful in Table [TUl 

6 The metric TNoI 

When we apply Grobner techniques to a set of equations (either by calculating 
a basis or a normal form) we are, in some sense, trying to simplify the set 
of equations. In a zero-dimensional ideal, as shown in the Gianni-Kalkbrener 
Theorem |Gia89IKal89| . a purely lexicographic Grobner basis has a very distinct, 
triangular structure. 

With this in mind we thought it may be of some use to consider the number 
of variables present in a certain problem and so defined the following quantity, 
TNoI, which stand for "Total Number of Indeterminates" : 

TNol(^^) = (9) 
where NoI(/) is the number of indeterminates present in a polynomial /. 



Table 9. IBH91I: effect of orderings =gc versus =gr 









=g/< 


ziR 






Time Cells 


Time Cells 


Intersection A 


c 


29426 


3763 


2470 


273 




R 






> 1000s 


? 


Intersection R 


c 


36262 


2795 


1482 


189 




R 






> 1000s 


? 


RcLndom A 


c 


17355 


1219 


570 


165 




R 






> 1000s 


? 


R,cLndom B 


c 


356670 


7119 


470 


141 




R 






> 1000s 


? 


Ellipse A* 


C 


262623 28557 


62496 14439 




R 






271726 29939 


Ellipse B* 


C 


> 1000s 


7 


> 1000s 


? 




R 






> 1000s 


? 


Solotareff A* 


C 


16014 


1751 


2025 


297 




R 






> 1000s 


? 


Solotareff B* 


C 


43439 


6091 


1647 


243 




R 






> 1000s 


? 


Collision A* 


C 


216028 


7895 


> 1000s 


? 




R 






> 1000s 


? 


Collision B* 


C 


> 1000s 


? 


> 1000s 


? 




R 






> 1000s 


? 



We note that =gr is definitely worse than =gc- 



Table 10. Spheres and Cylinders: <An. — degrees 





<zlR 

degrees Time Cells 


~g/<aii 
degrees Time Cells 


=g/^ /<a-r 
degrees Time Cells 


Si,S2,C 
&, S3, c 
S3, S4, c 


6 / 18 8654 1073 
6 / 19 189202 12097 
6 / 21 248340 11957 


5 / 9 905 267 
5 / 11 5911 1299 
5 / 15 8159 1359 


5 / 7 270 99 
5 / 10 499 213 
5 / 15 580 213 



'degrees' is td{An) / sotd{A„) . 



Table 11. [BH91) : degrees 







<ZiR 




=g/<aii 






degrees 


Time Cells 


degrees Time Cells 


Intersection A 


6/14 


29426 


3763 


17 / 50 2470 


273 


Intersection B 


6/14 


36262 


2795 


15 / 41 1482 


189 


Random A 


9/16 


17355 


1219 


19 / 68 570 


165 


Random B 


9/16 


356670 


7119 


19 / 73 470 


141 


Ellipse A* 


6/24 


262623 28557 


6 / 26 62496 14439 


Ellipse B* 


6 / 24 > 1000s 


? 


25 / 253 > 1000s 


? 


Solotareff A* 


10 / 25 


16014 


1751 


10 / 28 2025 


297 


Solotareff B* 


10 / 25 


43439 


6091 


21 / 69 1647 


243 


Collision A* 


6/23 


216028 


7895 


27 / 251 > 1000s 


? 


Collision B* 


6 / 23 > 1000s 


? 


36 / 875 > 1000s 


? 



6.1 TNoI data 



The results of calculating such a quantity are given in Table 8, Table 9 and 
Table 10, showing a promising correlation to whether our preconditioning (with 
compatible ordering) is beneficial or not. In particular we note the following 
points: 

— In every example where preconditioning reduces TNoI (15 cases) there is a 
significant reduction in timing (a decrease factor ranging from 4.20 to 757.26) 
and number of cells produced (a decrease factor ranging from 1.98 to 65.13). 

— When preconditioning increases TNoI (7 cases) then generally there is an 
increase in time (an increase factor ranging from 1.79 to 6.13) and the number 
of cells created (an increase factor ranging from 1.35 to 3.52) or the problem 
remains infeasible. There is one 'false positive' result f |CMMXY09[ Example 
2]) where there is an increase in TNoI but a slight improvement in the time 
(a decrease factor of 1.20) and cells produced (a decrease factor of 1.55). 

— TNoI alone does not measure the abstract difficulty of the calculations: In- 
tersection A has a higher TNoI than Ellipse A yet the latter takes 25 times 
longer and produces over 50 times as many cells. We have only shown how 
to use it to compare variants of the same problem. 

As mentioned above, calculating TNoI alone is not of a huge use, and even 
considering the difference or ratio does little to predict the degree of improve- 
ment to expect. However, if we take the logarithm of the ratio (equivalently the 
difference of the logarithms) of TNoI and compare to the time or number of cells 
we get some interesting results. 

Plotting these quantities against each other certainly suggested there was a 
positive correlation. Recall that the sample correlation coefficient is defined as 



and is a number between -1 and 1 that indicates how correlated data is. A 
sample coefficient of 1 indicates perfect positive correlation and a coefficient of 
-1 indicates perfect negative correlation. Although we are only working with 
a small bank of data (22 examples) and partially incomplete data (timings of 
> 1000s were replaced by 10000 seconds and unknown cell numbers were replaced 
by 100000 to allow for coefficient calculation) there were still promising results. 

Let S be the polynomial input, Vs its corresponding CAD, ts the time taken 
to calculate Vs and cs the number of cells in Vs- Let G be the Grobner basis 
calculated with respect to the compatible ordering and define Vq, to and cq in a 
similar fashion. With the data set we obtained the sample correlation coefficients 
were as follows: 

— comparing log(TNoI(5')) — log(TNoI(G)) with log(fs) — log(tG) gives a sam- 
ple coefficient r — 0.821 which indicates strong correlation (for our limited 
sample set). 



rx,Y = 



Y:=i{X.-X){Y^-Y) 



(10) 



— comparing log(TNol(S')) — log(TNol(G)) with log(cs) — log(cG) gives a sample 
coefficient r — 0.829 which again indicate a strong correlation (for our limited 
sample set). 

Of course correlation does not imply causation, especially with a relatively 
small data set, so let us look more deeply at what TNoI is measuring. 



Table 12. TNoI for Spheres 





<AH. 

TNoI Time Cells 


~g/<aii 
TNoI Time Cells 


=g/^ /<ar 
TNoI Time Cells 


51, S2, c 

52, S3, c 

53, S4, c 


8 8654 1073 
8 189202 12097 
8 248340 11957 


5 905 267 

6 5911 1299 

7 8159 1359 


4 270 99 

6 499 213 

7 580 213 



Table 13. TNoI for |BH91] 







<AR 






~g/<ar 




TNoI 


Time Cells 


TNoI 


Time Cells 


Intersection A 


8 


29426 


3763 


7 


2470 273 


Intersection B 


8 


36262 


2795 


7 


1482 189 


Random A 


9 


17355 


1219 


5 


570 165 


Random B 


9 


356670 


7119 


5 


471 141 


Ellipse A* 


7 


262623 28557 


6 


62496 14439 


Ellipse B* 


7 > 1000s 


? 


21 > 1000s ? 


Solotareff A* 


9 


16014 


1751 


8 


2025 297 


Solotareff B* 


9 


43439 


6091 


7 


1647 243 


Collision A* 


7 


216028 


7895 


18 > 1000s ? 


Collision B* 


7 > 1000s 


? 


22 > 1000s ? 



6.2 What is TNoI measuring? 

Consider what causes TNoI to decrease. Let 5' be a set of polynomials in variables 
Xi, . . . ,Xn ordered xi < X2 < ■ ■ ■ < Xn- The following are three possible reasons 
for a decrease in TNoI: 

1. The number of polynomials in a specific set of variables, {xi^, . . . ,Xi,}, is 
decreased. If Xk is the most important variable then reducing the number of 
these polynomials will simplify the decomposition in the (xi, . . . ,a;fc)-plane. 
This will simplify the overall CAD, reducing the number of cells produced 
and hence the time taken to calculate the decomposition. 

2. At least one variable is eliminated from a polynomial. If the variable Xk 
is eliminated from a polynomial p then the decomposition based around p 



Table 14. TNoI for |CMMXY09| 







<ZiR 




~G/<AYi. 






TNoI 


Time Cells 


TNoI Time 




Cells 


Cyclic-3 


9 


3136 


381 


6 20 + 245 = 


265 


21 


Cyclic-4 


16 > 1000s 


7 


6 64 + 5813 = 


5877 


621 


2 


7 


2249 


895 


14 22 + 1845 = 


1867 


579 


4 


6 


3225 


421 


11 24 + 19738 = 


19762 


1481 


6 


4 


363 


41 


5 20 + 918 = 


938 


89 


7 


8 


3667 


895 


22 26 + 6537 = 


6563 


1211 


8 


6 


3216 


365 


5 21 + 174 = 


195 


51 


13 


9 


14342 


4949 


4 18 + 220 = 


238 


81 


14 


11 


334860 27551 


9 21 + 971 = 


992 


423 



will be greatly simplified. This will again simplify the overall CAD, reducing 
the number of cells produced and hence the time taken to calculate the 
decomposition. 

3. A polynomial in a large number of variables, say fc, is replaced by j poly- 
nomials each with n.i variables such that ^ rt^ < fc. Intuitively this would 
increase the number of discriminants and resultants calculated, be it in the 
projection phase of <coi or in ^ac, but the results appear in lower levels of 
the projection tree, and this effect is more potent than the apparent increase 
in the number of discriminants and resultants. We have yet to build a good 
model of this, though. 

Obviously, in general, a combination of these factors will be the reason for the 
decrease in TNoI. Also, there may be opposing increases in TNoI, which pre- 
sumably explains why the 'false positive' of |CMMXY09( Example 2] shows an 
increase in TNoI but an improvement in the CAD efficiency. 

7 Conclusions 

— For both <coi and <aii and ^AC^ pre-conditioning the equations (where 
applicable) by means of a Grobner calculation is often well worth doing. 

— Grobner reduction of inequalities with respect to equalities has never, in our 
examples, made things worse. 

— A priori^ it can be quite difficult to see which combinations of Grobner base 
and Grobner reduction will be best, but the Grobner side is generally chea]:0. 

— We therefore have multiple equivalent formulations of a given problem. We 
have investigated the metrics of |DSS04) . but have concluded that, at the 
level of choice of formulation, TNoI is a better predictor. It does not help for 
predicting the best ordering of variables, for which }DSS04| or the Brown 
heuristic jBro04| are appropriate. Phisanbut |PhilH Chapter 8] found the 
Brown heuristic sufficiently good, and simpler to compute. 

* This is a significant change from |BH91] . who had examples where the Grobner cal- 
culations was much more expensive than the Cylindrical Algebraic Decomposition. 



— In Section 13.21 we saw how g = A / > could be reduced to g = A 
sprecond(/, g) > 0. In principle, given si = A • ■ • A = A / > 0, after 
computing a Grobner base G for the s^, we could attempt a more general 
reduction of / by G. Pure NormalForm reduction has proved useful (Tables 
[6l[7]), but we do not have enough good examples to measure the utility of a 
more general pseudoremainder-like reduction. 
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